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Abstract 



We present preliminary results of particle-in-cell (PIC) simulations of collisionless mag- 
netic reconnection for conditions that apply to magnetically dominated pair plasma jets 
such as those in gamma-ray bursts (GRB). We find similar behaviour to previous authors 
but with reduced reconnection efficiency This results because we include the full electro- 
magnetic field dynamically and allow electric field fluctuations to grow. Although weak, 
these fluctuations impede early x-point formation by periodically decelerating and even 
reversing the inflow. 
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1 Introduction 



Presently, few testable theoretical predictions or observational diagnostics exist 
to robustly constrain whether GRB outflows are initially dominated by particles 
or Poynting flux, and what fraction of the particles are baryons vs. pair plasma 
(|Drenkhahn and Spruitl 120021 : iKoers and GianniosLl2007h . Whatever the case close 
to the engine, at large distances observations imply that particles are accelerated to 
highly relativistic velocities and dominate the energy budget. Magnetic reconnec- 
tion is a plausible mechanism of particle acceleration in GRB jet plasma if the free 
energy is carried by magnetic fields. 
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Reconnection has mainly been studied in the lab, the solar corona, and particularly 
the Earth magnetosphere, which produced the 'Geospace Environmental Model- 
ing (GEM) reconnection challenge': Here a simplified reconnec tion problem wa s 



studied and compared using several MHD, hybrid and PIC codes (Bir n et all 12001). 
Much less attention has been paid -either analytical or numerical- to the fundamen- 
tal reconnection mechanisms in either the relativistic regime or for a pair plasma. 

Reconnection in a (relativistic) pair plasma is very different from an electron-ion 
plasma. For the latter, the difference in skin depth between the electrons and ions 
and the corresponding scale lengths of the current sheet and reconnection region 
play an important role; reconnection is facilitated by Whistler modes, which can be 
related to the Hall current in the non-ideal generalized Ohm's law. In a pair plasma, 
however, the contributions from the electron and pos itron equations of motion t o 



the Hall current cancel and Whistler waves don't exist (IBlackman and Field 1 1993). 
The non-ideal term in Ohm's law that allows field lines to reconnect is an effective 
enhanced resistivity related to the electron and positron pressure tensors and the 
tearing mode is the dominant instability. 

To improve our understanding of magnetic reconnection for conditions applicable 
to GRB outflows we are performing si mulations with a ful ly explicit and highly par- 



allelized relativistic PIC code, Osiris (|Fonseca et all |2002|) . By tuning the plasma 



parameters from non to highly relativistic values and comparing results for a pure 
pair plasma and a plasma with an increasing proton fraction, we can obtain differ- 
ent particle acceleration spectra to be used to construct observational diagnostics 
of the jet's nature. This complements the existing PIC simulation work in the GRB 
outflow context which ha s focused on sho ck acceleration via the Weibel instability 
in ion-electron plasmas ( Med ved evl |2007l) . or magnetic pulse acceleration without 
magnetic reconnection (lLiangu2005|) . Spectra from each of these classes of simu- 
lation should ultimately be compared. 




Fig. 1. fii-#2 (left) and v\-V2 (right) vector fields during reconnection around one x-point 
on each current sheet, with £3 as background (color), for two different simulations. 

A simple geometry for studying reconnection is the Harris current sheet, where 
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two regions of nearly uniform magnetic field of opposite sign are driven into a thin 
current sheet by their own magnetic pressure. If the plasma can deviate from ideal 
MHD, the field lines can slip and reconnect at one or more jc-points, after which the 
reconnected field lines are strongly bent. When those field lines straigten to min- 
imise their tension, particles in the outflow region are accelerated. Free magnetic 
energy is thus converted to particle energy in the process. Using the Harri s configu- 
ration , we can compare to the results i n the GEM reconnection challenge (|Birn et al . 



200 1[) and its extension to pair plasma dBessho and BhattacharjeeLl2005|) . 



The simulations presented here are for a pair plasma and have two current sheets 
to facilitate use of periodic boundary conditions on all walls. In each half -plane we 
have an initial magnetic field normalized as: B = ^ y ieC J°p j tanh(x2/a)f i, where 
a = di/2 is the thickness of the current sheet and dj = c/ co p is the skin depth in 
terms of the plasma frequency (o p = ^Ann^e 2 /m e = c£l c / (\/2va) and we choose 
a non-relativistic) Alfven velocity va = Bq/ ' \J%%m e n§ — 0.05c. In a steady state 
(dE/dt = 0) we have j = —2enoVACOsh~ 2 (x2/a) and n e = n p = nocoshr 2 {x,2/ a) 
with thermal velocity spread (v t h) XJjZ = va/V2, such that p = In^kT = Bo/Sn. 
A small uniform particle density = 0.2no is added to replenish the accelerated 
particles leaving the reconnection region and the total number of particles is 10 8 — 
10 9 . The spatial and time resolutions are Ax = Ay = 0.05J, and co p At = 0.035 with 
typical runs exceeding t = 10 3 /(D p . Most simulations use a L x x Ly = 512 x 512 
grid of dimension — 12.8J, < xi^ < 12. 8J,, but for some we double the length along 
the current-sheet to reduce the effect of periodic boundaries. To quickly force the 
system into the non-linear regime, sometimes we impose an initial perturbation, 
B x = 10- 4 cos^-cos^M v = 10- 4 sin^sin^. 



3 Results 



Figure 1 (left) shows the magnetic field at Q. c t = 64 when the reconnection has 
reached a quasi steady state, and the reconnection electric field is well developed. 
This simulation was started from a perturbed equilibrium magnetic field. The elec- 
tric fie ld at th e x-points is ~ ±0. 16Bqva , w hich is slightly lower than in lBirn et a3 
(1200 If) and lBessho and Bhattacharjed (|2005l) . When we start from numerical noise, 
the secondary tearing mode quickly develops, and the magnetic field fragments into 
a series of small x- and o-points that move under the magnetic pressure and finally 
merge into one dominant x-point on each current-sheet. The right panel shows the 
in- and outflow velocities (V1-V2) around the jc-points. One can see the particles 



flowing in with inflow velocity from \E Z 



and being accelerated along the 



x\ axis. In this figure the effect of the periodic boundary conditions is already be- 
coming visible with accelerated particles re-entering the box, but they're not yet 
influencing the reconnection region. 
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4 Discussion 



In all simulations of non-relativistic reconnection we've seen in the literature, the 

fluctuating electric field is either neglected by not evolving the displa cement current 

(com mon in MHD codes), or artificially eliminated at each time-step (|Bessho and Bhattacharjed . 



2005). Osiris does neither and evolves the full fields, including all electromagnetic 



wave modes. Doing so reveals weak plane-wave like fluctuations in £3 (due partly 
to numerical noise, and partly to particles being allowed to cross the current sheet 
that don't cancel out exactly). When this oscillating component has the 'wrong' 
sign, the flow of plasma into the current sheet is impeded, and initially this com- 
ponent is strong enough to reverse the bulk plasma motion away from the current 
sheet. Only during the wave-phase where £3 has the proper sign do the particles 
flow into the current sheet. Eventually, the x-points and a mean steady reconnec- 
tion electric field develop and fast non-linear reconnection is achieved. The final 
in- and outflow velocities reached in our simulations are somewhat slower than in 
previous work, and the evolution to a steady-state takes a bit longer. 



5 Future work 

Currently, we are analysing larger runs on the NERSC clusters using more par- 
ticles, higher resolution, and larger aspect ratio along the sheet. Furthermore, the 
implemented distribution functions in Osiris are only valid in the mildly relativistic 
regime, so we're implementing the proper relativistic Jiittner-Synge distribution. 
Then we can proceed to vary the Alfven velocity arbitrarily, adding a small fraction 
of baryons, and 3D generalizations where possible. Finally, we will calculate the 
radiation spectra. If the synchrotron cooling time is long compared to the reconnec- 
tion rate, this can be done post-procees. If the cooling time is fast, but the plasma 
is optically thin to the produced radiation, we must calculate the radiation self- 
consistently at each time-step. Differences in the spectra produced by reconne c- 



tion in a magnetically dominated wind (AC model iDrenkhahn and Spruitl ([2002)), 
a magnetic pulse model, or shock- acceleration in a matter dominated model, will 
provide a useful observational diagnostic to distinguish between the three. 
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